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NOMENCLATURE 



SYMBOL 


MEANING 


a 


Sound speed (ft/sec) 


A 


Area (ft^) 


BR 


Bypass ratio 


Cp 


Specific heat at constant pressure (ft-lb/slug-°R) 


FAR 


Fuel-air ratio 


g 


Acceleration of gravity (ft/sec^) 


H 


Enthalpy (ft-lb/slug) 


HVF 


Lower heating value of fuel (BTU/lbm) 


J 


Joules constant (778 BTU/ft-lb) 


M 


Mach number 


MR 


Mixture ratio 


MW 


Molecular weight 


m 


Mass flow rate (slug/sec) 


P 


Pressure (Ib/ft^) 


PP 


o 

Partial pressure (Ib/ft ), (atm) 


R 


Gas constant (ft-lb/slug-°R) 


S 


Entropy (ft-lb/slug-'^R) 


SFC 


Specific fuel consumption (Ib^^j/lb-hr) 


ST 


Specific thrust (Ib^-sec) 

Ibm 

m 


T 


Temperature (°R) 


TH 


Thrust (lb) 
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SYMBOL 


MEANING 


V 


Velocity (ft/sec) 


V 


Mean velocity (ft/sec) 


WGR 


Water to gas generator air ratio 


X 


Specific humidity 




Ratio of specific heats 




Efficiency 


? 


Density (slug/ft^) 


SUBSCRIPT 


MEANING 


A 


Air 


B 


Burner 


C 


Compressor 


D 


Gas generator diffusor 


f 


Saturated liquid 


fg 


Change by evaporation 


fu 


Fuel 


F 


Fan 


FD 


Fan diffusor 


g 


Saturated vapor 


i 


Refers to property at station (i) , i = 1,2... 


I 


Refers to state reached by isentropic process 


N 


Gas generator nozzle 


P 


Propulsive 


REFA 


Reference of air 
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SUBSCRIPT 



MEANING 



S Static 

T Total, 

W Water 

V Vapor 



Turbine 
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CHAPTER I 



INTRODUCTION 

With the present day emphasis on nuclear propulsion and its far 
reaching effects on a modem Navy, advances in other types of marine 
propulsion may appear less necessary and certainly less glamorous. 
However, the significance of small to medium displacement craft in a 
limited conflict has been forcibly brought home by the role such craft 
are playing in Vietnam and in other troubled and far-reaching corners 
of the globe. One of the major features of the limited conflict is the 
speed and mobility of participating units, generally with respect to 
land and air forces; however, naval forces must soon follow with faster 
vessels. With existing types of propulsion, speed increases must be 
paid for with dramatic weight increases often outweighing the speed 
increase. Consequently, one of the most recent fields of endeavor in 
marine propulsion has been to develop a system that supplies enough 
power to attain high speeds and yet is lightweight by comparison. This 
paper deals with one of the most recent proposals, the water-augmented 
turbofan engine . 

Gas turbines are currently being used to drive conventional screw 
propeller vessels and more exotic systems such as the water-jet driven 
hydrofoil (PCH-2) of the U. S. Navy. However, little serious thought 
has been given to using turbines for jet propulsive power since the high 
exit velocities inherent in gas turbines (even in turbofan engines) 
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produce characteristically low propulsive efficiency and are thus too 
impractical. Conventional thrust augmentation by water injection was 
considered; however, this generally involved injection of water into the 
compressor air flow or into hot turbine exhaust gases. In the case of 
the former, high purity water would be necessary because of the contact 
with vital moving parts. In both cases the high temperatures of the gas 
generator would vaporize most of the water, and the energy required to 
accomplish this would negate the thrust increases (1) . A mixed flow 
turbofan where fan duct and gas generator exhaust through the same 
nozzle would have the same deficiency since high temperatures in the 
combined nozzle would again vaporize most of the water. Davison and 
Sadowski (2) proposed a novel method of injecting sea water into the 
relatively cool flow of a fan duct where the fan and gas generator exhaust 
through separate nozzles. As a consequence, very little water is 
vaporized, and no moving parts are contacted, making the use of sea- 
water feasible. See Figure 1. 

The purpose of this paper is to optimize the design point parameters 
of bypass ratio, water to gas generator ratio, and those velocity, 
temperature, and pressure ratios that directly affect the performance of 
this system with respect to specific thrust and propulsi\e efficiency. In 
particular, it is desired to determine to what extent departures from 
thermodynamic and dynamic equilibrium in the mixing processes affect 
performance and what critical limits exist, if any. It would then have to 
be established experimentally that a mixing duct with two-phase flow in 
the actual conditions would satisfy these requirements. 
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WATER SPRAY NOZZLES 
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FIGURE 1 

WATER-AUGMENTED TURBOFAN ENGINE (ACTUAL) 



This is both necessary and important for two reasons. First, the 
work of Davison and Sadowski (2) , which constitutes the only major 
effort on the subject, deals with one particular engine and one fixed 
set of operational parameters. In other words, they have shown it to 
be feasible for one particular engine but have made no effort to find 
under what conditions such a system would be optimum, and such 
information is necessary for the design of a prototype. Second, an 
experimental mixing duct/fan nozzle section should be constructed to 
experimentally verify the theoretical findings of this paper, especially 
with regard to acceptable limits for velocity and temperature lags 
between water and air and also mixing duct lengths necessary to achieve 
them. 

The data obtained in the paper were obtained from a computer 
program designed to analyze the performance of the water-augmented 
turbofan engine while parametrically varying the critical design parameters 
so that graphical and tabular displays of thrust and efficiency variation 
could be obtained and optimum points determined. 

Briefly, the remainder of the thesis will present a description of 
the model used, develop the equations upon which the computer program 
is based, present methods of analyzing the data, and conclude with 
results, discussion, conclusions, and suggestions for further work to 
be done in this area. 
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CHAPTER II 



DESCRIPTION OF MODEL USED 

2 . 1 Physical Layout 

Figure 2 is a block-schematic diagram of the operating system 
depicted in Figure 1. Important points to consider are that the fan duct 
system is completely separate from the gas generator system with 
separate exhausts and that water injection occurs in the mixing duct 
entrance after the fan, thus contacting no moving parts. All stations 
are numbered in Figure 2; and when a numerical subscript is involved 
in an equation, the quantity thus subscripted is measured at that station 
in the engine. It should be noted that the figures do not show fuel tanks, 
w'ater injection pumping systems, or any other components contained in 
the actual engine. 

2.2 Assumptions 

Because of the large number of variables in an analysis of this 
type, many simplifying assumptions had to be made to limit computational 
time. The air was assumed to be behaving as a perfect gas with 
constant specific heats. As is common practice in rough design work 
in turbomachinery applications, Cp was assumed constant at one value 
for low temperatures and constant at a second value for high temperatures. 
Gamma was then computed from the appropriate value of Cp. For the 
purposes of this investigation, a plot of Cp versus temperature of the 
gas was made from data in reference (3). The results of this were two 
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FIGURE 2 

WATER-AUGMENTED TURBOFAN ENGINE ( BLOCK ) 



values for Cp for the two temperature conditions. For the combustion 
chamber, which is the transition between the two extremes, an arithmetic 
mean of the two values of Cp was used. This was justified on the basis 
that the plot of Cp versus temperature was roughly a straight line in the 
region of interest. Viscous effects were ignored with the exception of 
the influence of the water on the air in the fan duct. It was assumed 
that no heat was transferred other than in the combustion chamber, and 
no work was transferred other than in the turbine, compressor, and fan. 
Mixing duct area was assumed constant, and it was also assumed that 
the engine would have only converging nozzles. In other words, both 
gas generator and fan exit at atmospheric pressure; however, if this 
produced an exit Mach number greater than one, the values were 
recomputed to exit at a higher pressure with sonic velocity. The last 
major assumption made in the initial analysis was that the air entering 
the system was dry air with no relative humidity effects. For ease in 
calculations one slug of air was assumed to be entering the gas 
generator per unit time. As a consequence, bypass ratio was also the 
mass flow through the fan section per unit time, and an increase in 
bypass ratio keeping other parameters constant was an increase in fan 
section size holding gas generator size constant. 
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CHAPTER III 



ANALYSIS OF DRY TURBOFAN 



With zero water injection, the analysis of the system is developed 
from H“S diagram relationships and basic equations of engineering 
science. This chapter and the equations contained in it are arranged in 
sections similar to those contained in the computer program for the dry 
turbofan engine. 

3. 1 Gas Generator Diffuser 

Ambient total pressure is determined from 




( 1 ) 



so that defining diffuser efficiency as 




-Pr, - Pso 



( 2 ) 



it follows that 





(3) 



Since no work is done in the diffuser, 



1ti - Tro = Tso 1-^ 2. 



(4) 



and enthalpy is 




(5) 



where H 



REFA 



is the enthalpy of the air at T^^. 
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3.2 Compressor 



The ideal (isentropic) temperature at compressor exit for a specified 



pressure ratio 



Ptz 



is 




Tt» 



?T^ 

Pn 









/ 



( 6 ) 



Define compressor efficiency as the ratio of isentropic work input to 
actual work input; thus 






rz. 



c 



Pc 



Ttzi - Iti 






+- VAt \ 



T - 

I XL - : 'SO 



3.3 Combustor (Burner) 



Recall from Section 2.2 



^?c 



Cp„ = Icpc ^ 



Z.O 



fixed by maximum allowable turbine inlet temperature. 
Therefore 






Cfo Tt5 “ \tz Hti 



Writing the energy equation across the burner, 



(7) 

( 8 ) 



(9) 






m 






If it is assumed that injection temperature of the fuel is equal to air 
temperature at that point, and the energy needed to vaporize the fuel 
is included by choosing the lower heating value of the fuel, then 
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(10) 



FAR 



I 'T3 - I 

J - C?,,(t-,5-Tt2.| 



3.4 Turbine 

For the entire system work input must equal work output; therefore, 
may be calculated from 



m. 



and 



F\t4- ~ Ht3| “ V\ti Hti 






Htt ' V4t< 



Tt3 - 



( 11 ) 



( 12 ) 



Turbine efficiency is the ratio of actual work output to isentropic work 
output, so that 

Tt« = Tt5 - ~ M (13) 

'Hr 

An expansion whether isentropic or actual is to the same pressure; 
therefore , 



?T4 - - 






ITAI 

Tt3 






(14) 



3.5 Gas Generator Nozzle 

Since no work is done in the nozzle, 

Tta- - 



Knowing that the flow exits to ambient pressure fixes the pressure 



ratio 



2sjl 



Thus 



'l^-\ 




( 15 ) 
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Define a nozzle efficiency 



(16) 



Vi = 
In 



Tt4 - \ 



^ Tt4 - 



Equation (16) may be solved forTgg and exhaust velocity determined 
from 



Va 5 = - la's) 



( 17 ) 



Sound speed is 



Q-5 - 



Therefore exit Mach number 



K 



Vi\s 

(3.S 



(18) 



Note the comments in Section 2.2 concerning procedure when I>-1 



3.6 Fan Diffuser 

Specifying a separate fan diffuser efficiency , the equations 

for the fan diffuser correspond to equations (2) through (5) with subscript 
6 substituted for subscript 1 . 



3.7 Fan 



With a corresponding fan pressure ratio 



? 



T7 






and fan efficiency. 



P , the fan equations are identical to the compressor equations of 
Section 3.2 with the subscripts 7 and 6 substituted for the subscripts 
2 and 1 respectively. 



3.8 Fan Mixing Duct 

No work is done in the mixing; therefore. 
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(19) 



\T7 












?T7 



( 20 ) 



3 . 9 Fan Nozzle 



The equations of Section 3.5 are valid for the fan nozzle also with 
subscripts 9 and 8 substituted for 5 and 4 respectively and the proper 
fan nozzle efficiency specified. It should be noted at this point that the 
entropy at all stations is calculated from relationships similar to 

Tc '•n f) ( (■ 



St: - ^ T. 






~ 'RUrv 






^ s 



Re FA 



(21) 



where is the entropy at TgQ and Pqq 



3.10 Specific Thrust, Propulsive Efficiency, and Specific Fuel 
Con sumption 

Thrust is determined from the momentum flux equation 

jy\ = ^ rnpNK=\-Vo') + ^5(v^ss ' (22) 

Specific thrust based on the total mass flow of air entering the system 
is 

TH 

ST" " + (23) 

Propulsive efficiency is a measure of how well the system is converting 
the change in kinetic energy of the working medium into thrust power of 
the vehicle and can be shown to be 
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( 24 ) 




I I VA • Vq I 

TH • voj -V rn y( ■*" ITrAtr(V/^c\ - \/o| 

' ' 2: 2^ 



Specific fuel consumption is the number of pounds of fuel required to 



produce one pound of thrust continuously for one hour and can be shown 
to be 



SFC = FAR • ~ ^ ( 25 ) 

TH 
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CHAPTER IV 



FIRST WATER INJECTION ANALYSIS 

The injection of water in the mixing duct does not affect the results 
of Sections 3.1 through 3.7. The equations contained in these sections 
are still valid for the case of water injection. The new analysis is 
based on solving the problem in the mixing duct and fan nozzle and 
evaluating the new specific thrust and efficiency from the information. 

In this first water injection analysis it was assumed that none of 
the water in the mixing duct and fan nozzle sections was vaporized. 

The rationale behind this assumption is that, according to reference (2), 
much less than one per cent of the water is vaporized. 

4 . 1 Definitions 

The following terms will bo used throughout the chapters dealing 
with water injection: 

Bypass Ratio (BR) = vn f c 

Water to Gas Generator Ratio (WGR) = FOc 

Mixture Ratio (MR) = f 

The mean velocity at a station is defined after Elliot (4) as the 
velocity that the air and water would attain if they were allowed to 
reach equilibrium at that station with no loss in momentum. 



2.4 



thus , 

Mean Velocity (V) = \Jh ^ HR • Vw ( 27 ) 

4.2 Mixing Duct 

The momentum equation for the assumptions of the paper is 



rap 






- - Ad 

The area occupied by the air is m p ^ p ^ 

The area occupied by the water is \j 

Thus , the total flow area is A = 



(28) 



I 



HR 



(29) 



‘ PwV/w 

Combining equations (29) and (28) and utilizing the fact that the mixing 
duct is assumed to be of constant area, equation (28) can easily be 
integrated and becomes 



Vg, = y~f 

The energy equation is 



■ ”?'38 



\ -V 



HR 



f K Va Pw 



(30) 



6 



fr\ 



S7 



-V 



VA7 

z 



rrs. 



\AJ I 






1 \Ase> ^^2) 


1 VAvSlg, ^ 




1 z 



(31) 



Equations (30) and (31) have five unknowns: PS8, VA8, VW8, TS8, and 
TW8 . With three equations lacking for a solution, it was necessary to 
specificy ratios of unknowns to serve as equations. These three 
remaining equations were 
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TS^8Wy\ 



TW&/TS6 



(32) 



VR8WA 



VW£>/\yA3 



(33) 



PSRS7 



?SS/?S7 



(34) 



These admirably suited the eventual task of determining the effect 
of deviations from equilibrium in this mixing process on overall 
performance . 

4 . 3 Fan Nozzle 

The momentum equation cannot now be integrated to yield an equation 
of the form of equation (30) since variation of area and pressure in the 
nozzle is unknown. Several attempts were made to provide a second 
equation either from known facts or from assumptions concerning area 
variation in the nozzle with varying results. The problem was rendered 
academic when the underlying assumption of no water vaporization was 
found to yield erroneous results. (See Section 4.4.) 

4 . 4 Conclusion of First Attempt 

A continuity equation was introduced into the mixing duct 
calculations in order to reduce the number of variable ratios specified. 

It was introduced in the form of . 




A 



or more specifically 





( 35 ) 
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The problem was now overspecified with more equations than unknowns. 



When two of the three ratios in equations (32) -(34) were used with 
continuity, momentum, and energy as defined in equations (30), (31), 
and (35), ridiculous values were produced for the unknown variables. 
Careful analysis of the equations involved revealed that they were correct 
for the given assumptions indicating that the assumptions needed 
re-evaluation. The net result was that the assumption concerning no 
vaporization of the water was not correct, and, indeed, vaporization of 
the water was found to play a large role in determining the solution. 
Admittedly a very small percentage of water is vaporized, but for two 
reasons its effect is significant. First, the entalpy of water \apor is 
roughly one-thousand BTU's per Ibm. greater than that of liquid water. 
Second, the amount of water present with respect to air (MR) is a number 
ranging in value from roughly ten to four hundred. Thus the contribution 
in an energy equation of vaporized water had to be taken into account 
along with the resultant change in water mass flow. 
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CHAPTER V 



SECOND WATER INJECTION ANALYSIS 



The second analysis is similar to that of Davison and Sadowski (2) . 
However, it differs in that their work involves the assumptions of thermal 
and dynamic equilibrium in the mixing duct, dynamic equilibrium in the 
fan nozzle, and zero pressure change in the mixing duct while this 
paper takes into account possible variations in these quantities. 

5 . 1 Steam Table Properties 

It was necessary in the second analysis to calculate the properties 
of enthalpy of the liquid (Hf) , enthalpy change by vaporization (Hfg) , 
enthalpy of the saturated vapor (Hg) , entropy of the liquid (Sf) , entropy 
change by vaporization (Sfg) , entropy of the saturated vapor (Sg) , and 
partial pressure of the vapor (PPv) . These are all complicated functions 
of temperature. A least squares method was used to fit data points 
read in from the steam tables (5) in the range of interest, and it was 
found that a cubic approximation gave accuracy down to the last 
significant figure in the steam tables. The exception was partial pressure 
of the vapor, the actual equation of which was simple enough to be 
included in the program directly. Since 






(36) 



and 




( 37 ) 
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it was necessary only to compute the properties Hf, Hfg, Sf, and Sfg, 
thus determining Hg and Sg. The equations generated by the four least 
squares approximations and the PPv calculation are 
HF = -5.2244716x10^ + (1 . 1 71 9422) (Tw) - (3 . 1882153x10“^) (Tw)^ 



+ (1.9421466x10 ^)(Tw)^ (38) 

HFG = 1.440548x10^ - (1 . 1017804) (Tw) + (1 . 0967588x10“^) (Tw)^ 

-(7.4163358x10“^) (Tw)^ (39) 

SF = -1.6947235 + (5 . 325xl0“^) (Tw) - (4 . 7499752xl0”®) (Tw) ^ 

+ (1. 8842757x10“^) (Tw)^ (40) 

SFG = 8.812883 - (2 .442480x10“^) (Tw) + (2 . 8390143x10“^) (Tw)^ 

-(1. 2441932x10“®) (Tw)^ (41) 



Equations (38) through (41) are valid in the range of temperatures 
from 510 °R to 660 °R. Tw, the temperature of the liquid, is in degrees 
Rankine, and the units of output for the curve fit are BTU/lbm for enthalpy 
and BTU/lbm-°R for entropy. In the equations that follow it will be 
assumed that Hfg, for example, has already been converted into ft-lb/ 
slug. 

log 10 (PPv) = log 10 (Pc) - | (^^) 

where 

Pc = 218 . 167 atm . 

X = (647.27 - T) °K 
T = (Tw - 492)(9/5) + 273. 16 °K 
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A = 3.2437814 



B = 5.86826x10 



-3 



C = 1.1702379x10"® 



D = 2.1878462x10 



-3 



and the answer, PPv, will be in atmospheres. Again, in the equations 



5 . 2 Mixing Duct 

The specific humidity, X, of the mixture is defined as the ratio of 
the mass of water vapor to the mass of dry air in the system. From 
perfect gas equations for dry air and water vapor, it can be shown that 



It was assumed that vapor temperature was equal to liquid temperature. 
Taking into account vaporization, the energy equation is now 



The momentum equation (29) and its integrated result (30) are still 
valid for the case of vaporization. The problem is now one of three 
equations: (44), (43), and (30) with six unknowns: VW8, VA8, TW8, TS8, 
PS8, andX8. By again specifying ratios of unknowns, equations (32) 
through (34), a solution is possible. However, an explicit solution is 



that follow, it will be assumed that PPv has been converted into Ib/ft^. 




(43) 
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not possible since equations (43) and (44) contain involved functions 



of TW8 through the enthalpy and partial pressure terms. VA8 and VW8 
may be solved using equations (30) , (33), and (34) . The remaining 
equations may be manipulated in such a manner as to yield a function 
of TS8 only (called FUN) which is equal to zero. A Newton Rhapson 
iterative solution may then be applied 



and DFUN is the derivative of the function with respect to TS8, both 
evaluated at TS8j_j^. When the correct value of TS8 is reached, FUN will 
equal zero, and TW8 and X8 may be calculated from equations (32) and 
(43). 

5 . 3 Fan Nozzle 

Along the same reasoning as in Section 5.2, specific humidity at 
the fan nozzle exit is 



Note again vapor temperature was assumed equal to liquid temperature. 

With varying area, the momentum equation cannot be solved unless 
a geometric configuration is specified. Therefore, the assumption was 
made that the net change in entropy through the nozzle is zero. That is 



TS8. = TS8 



- FUN/DFUN 



(45) 





(46) 





'*J^VcK 



(47) 
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Expanding 



np) 



^?c 



\sq 



- 



?sq 

^S8 



+ 









— 




^w-| - rNp 






- O 


(48) 



Since exit pressure is atmospheric, the equations (48) and (46) have TS9, 
TW9, and X9 unknown. By specifying 



TSR9WA = TW9/TS9 



(49) 



and solving equations (46) and (48) for a function of TS9 equal to zero, 
the Newton-Rhapson iterative method again yields values for TS9, TW9, 
and X9 . The energy equation may be written 



^ <= j j 

■V )^e> [ I 



rn uo 

i 




Y.c, rn p jj 



X 

Vwq 

X 




(50) 



Let VR9WA = VW9/VA9 



(51) 



and solve equations (50) and (51) for VW9 and VA9 










* 'sqj + Xg, - Xq + 



MR 






XgiveawA^ 1 
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1 


1^ \ -V Xq 


HR - 


•1 


/ Z (52) 
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5.4 Specific Thrust and Propulsive Efficiency 



Thrust, specific thrust, and propulsive efficiency are defined as 
before, but account must be taken of the varying mass flow rates in the 
mixing sections. 

TM = Va 5 - r^«c Vo ^ M ^ P I - mpVo 

"h I ho wj P I -* rn Vs)"/ VvAj o | | 



ST ^ 



1 -V FAR] 



Vas 



- )(q H>r \ 
\ 4 I 









\ + -&R, 



Vo 



\ 4 



V\M( 



4 - 



F 



[ e>R|\4Xqi 

I' BR| 

\4FAR 

Mbr) (=asVa 




•l p - TH • Vt^ [t H * Vo] 



o 

m 



T ( Vas ~ vJ 
2. 




(55) 



It should be noted that ambient air velocity was assumed equal to ambient 
water velocity. In other words, wind and current were neglected. 
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CHAPTER VI 



RESULTS AND DISCUSSION 

6 . 1 Organization of Output 

The results obtained focus on two basic indicators of water-augmented 
turbofan performance. These are the ratio of specific thrust for the water- 
augmented engine to specific thrust for the dry engine and the ratio of 
propulsive efficiency for the wet engine to propulsive efficiency for the 
dry engine. Hereafter these two parameters will be referred to as thrust 
ratio and efficiency ratio. The final data concern two broad areas of 
interest, the effect of variations in general system parameters on thrust 
ratio and efficiency ratio holding mixing section equilibrium parameters 
fixed at 1.0, and the effect of variations in mixing section equilibrium 
parameters on thrust ratio and efficiency ratio holding all other system 
parameters constant. 

Subsequent figures show plots of these variations and also the 
effect of water injection velocity on system performance. 

6.2 General System Parameters vs. Thrust and Efficiency Ratios 

The parameters originally thought to have the most significant effects 
on thrust ratio (TR) and efficiency ratio ( *1 R) are water to gas generator 
air ratio (WGR), ship speed (VO), and fan pressure ratio (FPR). Performance 
is also improved by increasing the bypass ratio (BR) , but this parameter 
was not thoroughly investigated since its limit is a function of design 
considerations only. 
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Figure 3 presents thrust ratio versus WGR for varying VO and FPR. 
The more water injected, the greater is its influence on the air as 
evidenced by a decrease in fan nozzle exit velocity (VA9) . From the 
thrust equation , equation (63) , it can be seen that a point is reached 
with increasing WGR where the negative contribution of the momentum 
drag of the water outweighs the positive contributions to thrust. That is, 
the momentum flux into the system eventually increases more rapidly 
than the momentum flux out. Thus, in general, the curves increase, 
reach a maximum or peak value, and then fall off as more water is added. 
Note that for a higher value of VO, values of TR are much lower. This 
is because of the increase in inlet momentum drag of the air and water 
with respect to other thrust contributions. Also note that at the higher 
values of VO, the curves peak at a lower WGR. Again, this is because 
the inlet momentum drag of the water predominates over the increase in 
thrust much more rapidly. 

In comparing Figure 3 with the curves obtained by Davison and 
Sadowski (2) , close correlation is found in shape and numerical \alues 
with differences attributable to the fan pressure ratio. In this paper fan 
pressure ratio is defined as the ratio of total pressure aft of the fan to 
total pressure before the fan, thus making mixing duct static pressure a 
variable with ship speed, while Davison apparently defines fan pressure 
ratio as the ratio of static pressure aft of the fan to ambient pressure, 
since his mixing duct static pressure is constant for all ship speeds. 
Figure 3 shows the fan pressure ratio to be an extremely important design 
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parameter when water injection is considered. The values of thrust 
ratio obtained clearly indicate that excellent values of peak thrust 
augmentation could be obtained in the order of upwards of three times 
dry thrust to over one and one-half times dry thrust. 

Figure 4 presents efficiency ratio versus the same independent 
parameters with a constant fan pressure ratio. The shape of the curves 
obtained is nearly identical to those in Figure 3 with roughly similar 
numerical values. Efficiency ratios peak at slightly higher values of 
WGR than the corresponding thrust ratio indicating that efficiency ratio 
is slightly less sensitive to the increase in momentum drag of the water. 
No particular peak value predominated o\er the range of velocities for 
either thrust ratio or efficiency ratio. However, for a specified \elocity 
and fan pressure ratio, a nearly identical value of optimum WGR was 
found for both thrust and efficiency ratios. Significant was the fact that 
at higher velocities this value was relatively low indicating that high 
speed craft are equally if not more suitable for water injection that those 
operating at lower speeds. 

6.3 Equilibrium Parameters vs. Thrust and Efficiency Ratios 

For the analysis of variations in mixing duct and fan nozzle 
equilibrium parameters on thrust and efficiency ratios, craft speed was 
fixed at fifty knots, fan pressure ratio at 1.5, bypass ratio at 4, WGR 
at 50, and compressor pressure ratio at 13.8. The four equilibrium 
parameters chosen were the ratio of water to air temperature at mixing duct 
exit (TW8/TS8) and fan nozzle exit (TW9/TS9) , the ratio of water to air 
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velocity at mixing duct exit (VW8/VA8) , and fan nozzle exit (VW9/VA9). 
Four runs were made holding all but one of the equilibrium parameters 
fixed at 1 . 0 and varying the remaining one through its maximum possible 
range, (Figures 5 through 8). For a slightly more meaningful output the 
two temperature ratios were varied simultaneously with the velocity 
ratios held constant (Figure 9) and likewise the velocity ratios varied 
together with temperature ratios fixed (Figure 10). 

Figures 5 and 7 indicate that for all other equilibrium parameters 
fixed at unity, equilibrium at the mixing duct exit is either detrimental 
or has no effect on thrust and efficiency ratios. This result was at first 
thought erroneous but then rationalized as being correct (although 
unrealistic). Note that for these curves, equilibrium conditions are 
required at the exit of the fan nozzle. Thus by requiring less equilibrium 
in the mixing duct, the mixing is being forced to take place in the 
nozzle, where it occurs isentropically , resulting in higher efficiency 
and thrust ratios . 

The results obtained in Figures 6 and 8 are much more significant 
as they represent not only physically possible states but, in the case 
of Figure 6, physically probable ones. Figure 6 shows that thermal 
equilibrium at the nozzle exit is definitely not desirable. Peak thrust 
occurs for values of TVV9/TS9 between 1.0 and 1.05, which can amount 
to as much as twenty-five or thirty degrees temperature difference. At 
greater values of this temperature ratio, thrust falls off again although 
much more gradually, indicating that a wise operating condition would 
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be to the right of the maximum point. Efficiency ratio is not affected 
to nearly the extent of thrust ratio although it also increases with a 
divergence from equilibrium and in fact never does begin to peak. 

Figure 9 plots thrust and efficiency ratios versu^emperature ratio 
at nozzle exit varying concurrently with temperature ratio at mixing duct 
exit. In this case, the predominance of the nozzle exit temperature 
ratio is again felt as perfect thermal equilibrium is seen to be the worst 
possible condition with a peak va lue occurring shortly after the 
equilibrium condition. Figure 10 shows the same type of plot with both 
velocity ratios varying and both temperature ratios held constant. It 
reflects the fact that the velocity ratio in the mixing duct makes little 
or no difference on thrust holding the other ratios fixed in that it is 
almost identical to Figure 8. Again departure from equilibrium means 
a decreasing thrust ratio and efficiency ratio. 

6.4 Other Results of Analysis 

One of the more significant results of the paper was discovered 
because of initial failures to assume realistic data. Initially a constant 
water injection velocity was specified regardless of craft speed. It 
was felt that this could be achieved with proper spray nozzle design and 
a pumping system to provide the necessary pressure differential when 
the speed of the vessel through the water did not create enough head to 
inject the water at the desired velocity. This constant injection velocity 
specification gave absurdly high values of thrust ratio not believed to 
be possible. With the new assumption that the injection velocity is a 
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fraction of the vessel speed (obviating the use of pumps) data were 
obtained which closely correlate with previous work (2) . However, the 
fact remained that, for a higher injection velocity, values of specific 
thrust were extremely high. Several TR versus WGR curves for varying 
injection velocity are presented in Figure 11 with equilibrium parameters, 
bypass ratio, and fan pressure ratio fixed. It was then postulated that 
if a pumping system were incorporated to provide the higher injection 
velocities, such a system could be used in emergencies and/or for 
important missions when a burst of additional thrust was required. 
Naturally, the power required to run such a pumping system would reduce 
the propulsive efficiency, but this would be well offset by the otherwise 
unobtainable increase in thrust. Note that this would not require more 
water injection but only higher injection velocities. In effect, it would 
be a type of sea-borne "afterburner." Such a device in intelligence 
ships such as the USS PUEBLO would be of unquestioned value for use 
in emergencies. Perhaps retractable hydrofoils could be incorporated to 
be used only when on the high thrust mode of operation. 



45 



14.0 



12.0 




100 200 300 400 500 

WGR 

FIGURE 11 

THRUST RATIO VS. WATER TO GAS GENERATOR RATIO 

46 



CHAPTER VII 



CONCLUSIONS AND SUGGESTIONS FOR FURTHER STUDY 
7.1 Conclusions 

The proposed concept of water injection appears to be a feasible 
means of increasing thrust and propulsive efficiency. Peak values of 
thrust ratio and efficiency ratio occur when momentum drag of the 
injected water predominates over the thrust due to additional mass flow. 
These peak values occur at lower values of water to gas generator ratio 
for increasing vessel speed because the momentum drag of the water is 
proportional to vessel speed. 

Peak values for thrust ratios occur slightly before peak values for 
efficiency ratios since the efficiency is less sensitive to inlet momentum 
drag of the water. 

Static pressure in the mixing duct plays an important role in the 
slope of the curves and thus determines the effectiveness of the water 
addition scheme at various ship speeds. 

The equilibrium parameter in the mixing process which appears to 
have significance in affecting system performance is the ratio of water 
temperature to air temperature at the fan nozzle exit. Peak performance 
with regard to thrust ratio is obtained when TW9/TS9 is in the range 
of 1.02 to 1.05. Considerably better values of thrust ratio occur on 
the side away from equilibrium as opposed to closer to equilibrium. 
Efficiency ratio continues to improve as TW9/TS9 moves further away 
from equilibrium. 
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Normally, injection velocity would be a fraction of the vessel 
speed determined by how well the ram scoop inlet is designed. If a 
pumping system is installed such that the water might be injected at a 
higher velocity, short periods of high thrust could be experienced with 
a corresponding decrease in propulsive efficiency, such as in an 
aircraft afterburner. 

7.2 Suggestions for Further Study 

There is a great deal of further work required before an experimental 
setup can be constructed', much of this work can be accomplished with 
this program and modifications to it. 

The ambient conditions were fixed throughout this analysis and an 
investigation into the variations of system performance for different sea 
water temperatures and air temperatures and pressures could be made. 
Since vaporization of the water in the mixing sections was of such 
significance, the effect of relative humidity in the air should be 
investigated. It should also be possible to choose a friction factor for 
the mixing duct and through the use of the Fanno flow equations 
determine the pressure drop in the mixing duct rather than specify it as 
was done in this analysis. The question of water injection velocity and 
a pumping system "afterburner" should be investigated with regard to the 
decreases in efficiency caused by the additional power required to run 
the pump. 

Last and most important, it remains to design and construct an 
experimental fan section that can investigate the effects of (a) the 
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equilibrium ratios, (b) static pressure in the mixing duct, and (c) injection 
velocities of the air and water on thrust and efficiency and to correlate 
the experimental results with the theoretical values obtained in this 
program and in this analysis. 
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APPENDIX I 



COMPUTER PROGRAM NOMENCLATURE 



Al.l Input Parameters 



Card 1: 



ETAD Gas generator diffuser efficiency 

ETAC Compressor efficiency 

ETAB Burner efficiency 

ETAT Turbine efficiency 

ETAN Gas generator nozzle efficiency 

ETAFD Fan duct diffuser efficiency 

ETAF Fan efficiency 

Card 2: 

DSO Ambient density (slug/ft^) 

PSO Ambient pressure (Ib/ft^) 

TSO Ambient temperature (degrees Rankine) 

HVF Heating value of fuel (BTU/lbm) 

TREFA Reference air temperature (492 deg. R) 



Card 3: 

AMN7 Air injection Mach number 
TWO Sea water temperature (deg. R) 

PTR32 Pressure ratio across burner 



Card 4: 



NVO 


Number of 


NCPR 


Number of 


NFPR 


Number of 


NTT 


Number of 


NBR 


Number of 


NWGR 


Number of 


NMDPR 


Number of 


NMDVR 


Number of 


NMDTR 


Number of 


NFNVR 


Number of 


NFNTR 


Number of 



vessel speeds 
compressor pressure ratios 
fan pressure ratios 
turbine inlet temperatures 
bypass ratios 

water to gas generator ratios 
mixing duct pressure ratios 
mixing duct velocity ratios 
mixing duct temperature ratios 
fan nozzle velocity ratios 
fan nozzle temperature ratios 
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Card 5: 



GC 


Cool temperature value of ratio of specific heats 


GH 
Card 6: 


Hot temperature value of ratio of specific heats 


VO (I) 
Card 7: 


Values of vessel speed (ft/sec) , I = l,NVO 


PTR21(J) 
Card 8: 


Values of compressor pressure ratio (Ib/ft^) , J = 1,NCPR 


PTR76(M) 
Card 9: 


Values of fan pressure ratio (Ib/ft^) , M = 1,NFPR 


TT3(K) 
Card 10: 


Values of turbine inlet temperature (deg R) , K = 1,NTT 


BR(L) 
Card 11: 


Values of bypass ratio, L = 1,NBR 


WGR(N) 
Card 12: 


Values of water to gas generator ratio, N = 1 ,NWGR 


PSR87(IJ) 
Card 13: 


Values of mixing duct pressure ratio, IJ = 1,NMDPR 


VR8WA(IK) 
Card 14: 


Values of mixing duct velocity ratio, IK = 1,NMDVR 


TSR8WA(IL) 
Card 15: 


Values of mixing duct temperature ratio, IL = 1,NMDTR 


VR9WA(IM) 
Card 16: 


Values of fan nozzle velocity ratio, IM = 1,NFNVR 


TSR9WA(IN) 


Values of fan nozzle temperature ratio, IN = 1,NFNTR 
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AI.2 Output parameters 



VA8 

VW8 

VA9 

VW9 

TS8 

TW8 

TS9 

TW9 

STHRUS 

AMN5 

AMN9 

ETAPRO 

SFC 

FAR 

X9 



Air velocity at mixing duct exit (ft/sec) 

Water velocity at mixing duct exit (ft/sec) 

Air velocity at fan nozzle exit (ft/sec) 

Water velocity at fan nozzle exit (ft/sec) 

Air static temperature at mixing duct exit (°R) 

Water static temperature at mixing duct exit (°R) 

Air static temperature at fan nozzle exit (°R) 

Water static temperature at fan nozzle exit (°R) 

Specific thrust (lb thrust/lbm/sec) 

Gas generator exit Mach number 
Fan nozzle exit Mach number 
Propulsive efficiency 

Specific fuel consumption (lb fuel/lb thrust-hr) 

Fuel-air ratio (lb fuel/lb air) 

Specific humidity at fan nozzle exit (lb vapor/lb dry air) 
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APPENDIX II 

COMPUTER PROGRAM AND SAMPLE OUTPUT 



rMPMAS C. KNLOS'^N 
,>)fT TURHOFAN 

sample data rum 

?.-■ \ 1 ARCHt 



real '^R 

dimension V0(ft).PTR?l(i:),PTR76(l''),TTS(S),HP(l'),WGRn'-»,PSPA7(i' 

- ) , \/R3WA( IJ ) , TSRRWA ( m , VRPi^AI 1 : ) , TSRDWAI 1 D ) 

STEAM TABLES FUNCTIONS AND DERIVATIVES 

HF(i<X)= -522.4A716 + 1 , 1 7 1 942 2 X - . ? 7 i rs?i 53*xx* XX ♦ .I'-'-'v-iq 

«^2 1A6F* X 

HFG( XX ) =I44-' . - 1,i:i78D4*XX 1 0 RS 7 SP * X V# X X - , E 74 ^ ^,3 

••:=358*XX*«3 

SF ( XX » =-l . 6RA7235 + j572sr'“^ I’i'XX - , r'C’ A '"'474Q0 7S 2* X X*X X + .''D-''-, 

■p D ' 0 I 8 B R2 7 S 7T-- X X #>-5 3 

Sr G( XX I =(f . 81 ?BB7 - . ' 24A24B0 I>p XX *■ . 1 7 R 3pr | 4 3 * X X * X X - , a '' 'V' " ^ ^ 
>!c2^41D 32’PX XX'>:-3 

0HF( XX) =I . 171Q422 - . )F, <7643 ; 6 >XXX + . 3 3 ' a T( R q? 6 39 8 X«X X« X X 

OhFG( XX )=- 1. If I 7804 + .V. 0 2I93SI 76*XX - , 3 ' f a -- 2 ? 240Ar a 74»rxx*x X 

OSH XX)=,OD 8 ^?SnOa) - .A!' ,,t^'D49P98'"4^XX <- . 3E 0 3 C D Ar q 6 6 7 8 2 ^ 1 * X X * X X 

DSFG( XX )=- . 3744248^1 ♦ . rOD'- 567R3286Xtxx - . 3< : ) f 3 ^ 3 7 7 ? 5 796 «X X* x v 

READ! 5, 1 ) ETAD, ETAC tETAB .ETAT tFTAN.RTAFO, FTAF 

1 FORMAT! 7F5. 2) 

RE AD( 8, 2)nS0,PSn,Tsn,HVF , TRFFA 

2 FnRMAT(F10.6,4Fl ). I ) 

READ! 5 , 3) AiMN7, TWC, PTP3 2 

^ FORMAT ( 3F8. 2 ) 

RE AD! 5, A ) NVOtNCRR t ^'FPO ,NTT tNPR , NWGR .NMUPR , NMDVR , NMDT R , NFNVR , NFNTP 
4 rORMAT! 1113) 

REAO(5,8)GCtGH 
8 FORMAT! 2F9.fc) 

Rf AD! 8 t 8 ) ! vn( 1 ) . I = I ,NVO) 

RF AD! 8, 6) ! PTF 21! J) t J = 1 tNCPR) 

RCAD!8,6)!PTR76!m),m=1,MFPR) 

R F AO ! 8 , 6 ) ! T T 3 ! XI t R - I t NT T ) 

RF AD! 8, 8) ! HR !l ) , L = 1 . NBR) 

RE AO ! 8 , 6) I WGP I N I ,N=I ,MNGR » 

READ! 5, 8) I PSE87! IJ ) , IJ = I , NmdpR ) 

READ! 8,8) ! VRPWA! IK ) , 1K = I, NMDVR ) 

REA0!8,6)ITSP8li,A!IL),IL = lt NMDTR ) 

READ!8,8)|VR9WA!1M),IM=1,nFNVP) 

RFAD!8,6) ! TSP9WA! IN) . IN=1 .NFMTP) 

8 PORMAT! RPP.2) 

R=I71 7. 81 
G=32.2 
CJ=778.26 
CPH=GH*P / ! GH- 1.3) 

CPC = GC“R / ! GC- 1. 7 ) 

TW7=TW0 

HRFFA= CPC*!TSn - TP^fA) 
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non onn ooo non ono 



IV G LEVtL 1| ^OD 0 



MAIN 



DATE = 6a0fll 



1 4/^7/S7 



HREFW=HE( TwaiOGJOG 
SREFA = r>,0 
WRITe(6,50I 
SO FOR AT ( • I ' I 
WPITE(6,7) 

7 FORMAT! //////T40, • TURBOFAN ENGINE CYCLE PROGRAM') 

WR1TE(6,8) 

B F0RMAT(//T4A,'WITH WATER INJECTION') 

WP ITE( 6,9) 

9 FORMAT! //T47, ' ALT I TUOE = 0 FT') 

WRITEI6, 10) 

10 FORMAT! ////T3, 'AMBIENT PROPER T I ES ', T35 ,' R ^FER E NC E PROPER T I ES ', T74, 
♦'EFFICIENCIES' ) 

WRITE!6,ll)OSn,HRFFA,FTAn,PSO,SREFA,FTAC,TSn,HREFW,FTAR,ETAT,ETAN, 
♦ETAFO, ETAF 

11 FORMAT! //T3, 'DSn = ' , F9. 6 , T2 1 , ' SL UG/CD. FT ' , T3 5 , • HR FF A= ',FH.1,T5 

♦6, 'FT-LB/SLUG' ,T74, 'FTAO = ' , F 6, 2 , / , T 3 , ' P SO = ' , F9, 3 , T21 , 'LB/SO, F 
♦T' ,T35, 'SRFFA= ' ,F 1 1 . 1 ,TS6, 'FT-LB/SLUG-OEG R',T74,'ETAC = ',F6,?,/ 
♦,T3,'TSn = ' ,F9,4,T21 , 'DEG R ',T35,'HREFW = ' , F 1 1 , 1 , T 56, ' FT-L B/SLUG 
*-OEG R',T74,'ETAR = ' , F6, 2 , / , T74 , ' FT AT = ' , F6 , 2 , / , T74 , • ET AN = ' , F6 
♦.2, /,T74, 'ETAFO = ' ,F5.2, / ,T74, 'ETAF = • , F6. 2 , / , T3 , ' WATER INJECTIO 

♦N SPECIFICATIONS') 

WRITE!6,12)AMN7,TW7 

12 FORMAT! //T3, 'AMN7 = ' , F 7 , 2 , / , T3 , ' VWO = VO KNOTS 

♦ ',/,T3,'TW7 = ’ ,F7.2,T19, 'DEG R*) 

no 3f) 1 = 1, Nvr 
CV=VO! I )♦! .6P8944 
VWO=OV 

AO=SORT(GC*R*TSOl 
AMNn=nv/AO 
vw7=o.B ♦ vwr 

GAS GENERATOR DIFFUSFR 

PTO = PSr*! !1.0+!!GC-l.C)/2.O)*AMN0*AMNP)>f-*!GC/!GC-l.r) ) ) 

PT1=ETA0 *!PTO-PSO) ♦ PSO 

TT1 = TS0*! 1 .f + ! !GC-1.0) /2. 1 ) ♦AMNO + AMNO ) 

ST1=CPC^AL0G! TTl /TSO)-R#ALOG! PT 1/PS0)+SRFFA 
HT1=CPC^! TTl-TSO) +MRFF A 

FAN DIFFUSER 

PT6=ETAFD=i'! PTO-PSO) PSO 

TT6 = TS0*! I.C + ! ( GC- 1 . " ) /2 . "^ ) ♦ AmnO* A mnO ) 

ST6=CPC*AL0G! TT6 /TSO) -R^ALOG! PT6/PS0) +SREF A 
HT6=CPC^! TT6-TS0)+HREFA 

GAS GENERATOR COMPRESSOR 

00 35 J=1 , NCPR 

TT2I=TT1*PTP21 ! J )♦♦! !GC-1.C ) /GO 
0HT21 I=CPC«! TT2I-TT1 ) 

HT2=DHT21 I /ETAC ♦ HTl 
TT2= ! HT2-HRFFA ) /CPC + TSO 
PT2=PTR2H J)*PT1 

ST2 = CPC.*AL0G! TT2/TS0) -R^ALOG! PT2/PS0) + SRFFA 
GAS GENFRATOP BURNER 
00 34 K=1,NTT 

PT3=PTR32*PT2 
CPB=(CPH 4- CPO/2.0 
HT3=CPR*! TT3!K )-TT2) + HT2 

ST3 = CPB«ALOG!TT3!K ) /TT2) - R*=ALOG! PT3/PT2 ) * SREFA 
FAR=1.0/( ( ETAB*HVF*G*CJ)/!CPH*!TT3!K)-TT2) )-1.0) 

FAN 

00 33 M=1 ,NFPR 

TT7I=TT6*! PTR76! «)*♦! !GC-1.0) /GO ) 

0HT76I=CPC*! TT7I-TT6) 

HT7=DHT76I /ETAF + HT6 
TT7=!HT7-MREFA ) /CPC 4- TSO 



noo ooo ooo ooo 



IV G I f^VEL !♦ MOD 0 



MAIN 






68081 



16 / 37/57 



PT7=PTR76( y) ^PT6 

S T7 = C PC « AL CG ( TT7/T SCI) -R^ALHG I PT7/PSn» + *nRFFA 
WPITF(6, 13) 

1 3 FORMAT ( / // // /# T3 f + + + + 

WV7=0. P*VC( I ) 

WRITE(6.14»VC( I > .PTR2K J) , PTR 76( M ) , TT 3 ( K ) , WV7 

14 FDRMATC ////Tin, • VO (KNOTS) IS • , F5 . 1 , / , T 10 , • C OMPRE S SOR PRESSUPF RA 
i^TIO IS ' ,F5. 1,/, TK, ‘FAN PRESSURE RATIO IS • , F5 . 1 , / , T 10 , • TIIR R I NF I 
♦ NLET TEMPERATURE (DEO. R) IS • , F 7. 1 , / , T 1 0 , ' WAT ER INJECTION VFI OC I T 
*y (KNOTS) IS • ,F5. 1, /////) 

GAS GENERATOR TURBINE 

DU 32 L=1,NBR 
BY=RR ( L ) 

CMF=l.n 

fmf=by 

TMF = CMF *■ CMF*FAR 

0HT43=(CMF*(HT2-HT1)+FMF*(HT7-HT6) ) /TMF 

HT4=HT3-DHT43 

0HT43I=0HT43/ETAT 

HT4I=HT3-0HT43I 

TT4I=TT3(K I-(HT3-HT4I ) /CPH 

TT4=TT3(K ) -( HT3-HT4) /CPH 

PT4I=pT3*((TT4I/TT3(K) ) ♦* ( GH/ ( GH- 1 . ) ) ) 

PT4=PT4I 

IF( ( PT4-PSD ,LT. ^.0) GO TO 24 

ST4=CPH*AL0G(TT4/TT?(K ) )-R*ALOG( PT4/PT3) + ST3 

GAS GENERATOR NGZ7LE 

PS5=Psn 

PS5I=PS0 

TT5I=TT4#( ( PS5I/PT4)»*( (GH-1.0)/GH) ) 

TS5=TT4 - ETAN*( TT4-TT5I ) 

TT5=TT4 

ST5=CPH*AL0G(TS5/TT3(K ) )-R*ALOG( PS5/PT3) ♦ ST3 
VA5=S0RT(2.C *CPH*( TT5-TS5) ) 

AS=SORT(GH*R*TSS) 

AMN5=VA5/A5 

STAnO=0.0 

IF( AMN5.LE.1.C) GO TO 15 
AMN5= 1. 0 

V A 5= SORT! ( 2.C*CPH>KTT5 ) /( 1 . 0‘-2 . 0*C Ph / ( GH*R ) ) » 

A5=VA5 

TS5=A5*A5/ ( GF*R ) 

TT5I =( FTAN*1 T4-TT4+TS5 ) /ETAN 
PS5=PT4*( TT5I/TT4)*«(GH/(GH-1,C ) ) 

ST5 = CPH*AL0G( TS5/TT3(K ) )-o*ALOG( PS5/PT^) ST3 
0S5=PS5/(R*TS5) 

STADO=( 1 .0+FAR ) /( ( 1 . G+ BY ) *DS5*V A5 ) ♦ ( P S5-PS0 ) 

FAN MIXING nOCT 

15 DO 31 N=1,NWGR 
MP=wGR ( N I /BY 
WMF=WGR (N)«CMF 

TS7=TT7/( l.t ♦(GC-l.n)/2.':*AMN7*AMN7) 

HS7= ( TS7-TS0)>KCPC + HRFFA 
A7=SDRT(GC*R*TS7) 

VA7=AMN7*A7 

VM7=( VA74-MRx=VW7) / ( 1.0 + MR ) 

PS 7 = PT7 / ( 1 .r-K GC-1 .f ) /2.'^t'AMN7«AMN7 )**(GC/(GC-1.0) ) 

DA7=PS7/(R*TS7) 

miking duct temperature iteration 

on 30 ij=i,NMnpR 
PSfl=PSRR7( I J ) Y0S7 
DW7=1.937B 
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AN IV G LEVtL 1, MOD C 



MAIN 



DATE 



6 80fll 



14/T7/‘^^ 



VM8=VM7 * ( PS7-PS8) * ( 1 . /( l. + MR n *( 1 . / (nA7’«'VA7) +MR/ ( 0W7*VW7) I 

DO ?9 IK = 1,N'M0VR 

VAH=VM 8 /( ( l.+MR*VR 8 WA( IK) I /( 1 , +MR» ) 

VW 8 =VR 8 WA( IK )*VA 8 
TS8=530. 

ZA= 3.243781A 
ZP= . 00986876 
ZC= .0OtCO0C11702379 
ZD^ ..■'02187RA62 
on 28 IL=l.NMnTR 

16 T)» 8 =TSRBV<A( IL)*TS 8 

ZF= CPC«( TT7-TS8 ) ♦ MR #( MF ( TW7 ) -HE ( TW 8 ) ) *C J*G ♦ MR*VW7*VW7/2. 

* - MR*vW 8 *V^ 8 /^ . - VA8*VA8/2. 

X=(37-4.11 - ( TSRflWAI IL)*TS 8 - 4Q2.)*5./9.) 

TK = (273.16 + ( TSR 8 WAI ILK^TSR - -^9?. ) #5. /9. ) 

PPV 8 = 1C . v*(ALQGin(?18. 167) - X /TK *( ( Z A + Z R *X+ Z C * X*>!^ 3 ) / ( 1 . 9 + 7 Q* v ) ) 
>S)*IA, 6959*144, 

Z =PPVR/ ( PSP-PPV 8 ) >^'^. 627 
PFIJN 8 = ZE/( HFG( TW 8 ) *CJ*G) - Z 
nZE=-CPC - mr*DHF( TW 8 ) *CJ*G 
0T = 5. /9,*TSRPWA( ID 
nx=- 5 , /Q.*TSP 8 kA( IL ) 

Z'^ = ((l. ♦zn*x ) *( ZR*OX ♦ 3, *ZC*nx*X**2 ) - (ZA ♦ ZP*X 4 - 7C*X*» 3 ) ♦( 70*nx 
* ) ) /( l.+70*X)**2 
ZN=(TK*()X - X*DT)/TK**2 

DPPV 8 = PPV 8 *(X/TK*ZM + ( ZA^ZR*X-t-ZC*X**3) /( l. + ZD*X)*7N) 

DZ=( (PS 8 -PPV 8 )* 0 PPV 8 + ppvH*nPPV8)/(PSP-PPV8)**2*^,67? 

DPFUN 8 = ( HFG( TW 8 ) *C J*G*DZE-7E*DHFG( T WR ) *C )#G) / ( HFGI TW 8 ) *CJ*G ) **7-D7 
TS 8 =TS 8 - PFIN 8 ZDPFUNS 
IF( PFIJNB.LE.C. lE-03) GH TO 17 
GO Tn 16 

17 TWR=T5P8WA( IL )*TSR 

X=(374,ll - (TSR3WA( IL)*TS 8 - 492.)*5./9.) 

TK=(273,16 + ( TSR 8 WA ( I L ) *TSR - 497,)*5./9.) 

PPV 8 = ir,«*(ALnGK( 218. 167) - X /TK ( 7 A+ z B*x + 7C*X ** 3 ) / II . 0 + ZD*X ) ) 

*)*14. 6959*144. 

XR=PPV8/(PS8-PPV8)*^.67? 

C 

C FAN NTZZLF 

C 

PS9-PSn 

TS9=520. 

no 27 IN= 1 ,NFNTH 
IR TlvO=TSR9WA ( I N ) *TS9 

70=(CPC*ALCG(TS8/TS9) - R * AL OG ( P SR/ ^59 ) + X R* S RG ( T WR ) *C J*G ♦ «»* 

* (SF( TWR)-SF( TW9) )*CJx<G) 

X=(374.11 - ( TSR9WAI 1N)*TS9 - 46 2 . ) ♦ 5 . / 9 . ) 

TK=(273.16 ♦ ( TSR9WAI I N»*TS9 - 492.)*5./9.) 

PPV9= ir . **(ALnG10I 218. 167) - X /TK *( ( Z A+ Z R*X 4 - Z C*X** 3 ) / ( 1 . f + 7 n* x ) ) 

* )*14.6959* 144. 

ZR = PPV9/( PS9-PPV9)*.622/TSR9WA( IN) 

PFUN9=Z0/ ( SFG( TW9 ) *CJ*G) - ZR 
nZ0=-CPC/TS9 - MR*nSF( TW9) *CJ*G*TSPPWA( IN) 
nr=5./9.*TSR9WA( IN) 

DX=-5./9.*TSP9WA(IN) 

ZM = ((1. +ZD*X) *( 7R*0X * 3. *7C*nX*X**2 ) - (7A -i-Z R*X ♦ 7C*X** 3 )*( 7 D*DX 
>^ ) ) / ( l.+ZD*X)**2 
ZN=(TK*0X - X*0T)/TK**2 

DPPV9= PP V9* ( X/TK* 7 M A- ( ZAfZR*X+ZC*X**3) /( l. + Zn*X)*ZN) 

DZ = ( (PS9-PPV9)*0PPV9 A- PPV9*DPPV9)/(PS9-PDV9)**?*0.622/TSR9WA(1M) 
DPFUN9=( ( SFGI TVi9 )*DZO-ZQ*DSFG( TW9)*TSR9WA( IN) ) *C, J*G) /( SFG( TW9) *0 J* 
*G )**2 - nz 

TS9=TS9 - PFLN9/DPFIJN9 

IF{ PFUN9.LE.''. lE-03) GO TO 19 

GC TO IR 

IQ TW9=TSR9WA( IN)*TS9 

X=(374.11 - ( TSR9WAII N)*TS9 - 49?.)*5./9.) 

TK=(273.16 ♦ ( TSR9WAI I N) *TS9 - 492.)*5./9.) 

PPV9= 1 ( .**( AL0G1CI218. 167) - X /TK *( ( Z A+ Z 8 *x* ZC*X**3 ) / ( 1 . 0+ Z 0*x ) ) 
’S')*14, 6959*144. 

xq=( PPV9/ ( PS9-PPV9 )*.622)/TSP9WA(IN) 
on 26 IM=1,NFNVP 

VA9=SQRT( (CPC*(TS8-TS9 ) ♦ X 8 *HFG ( TWR ) *C J*G * MR*HF | T W 8 ) *C J*G - X9* 
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ooo 



>1 IV G LFVEL 1, C 



MAIN 



DATE = 680.RI 



14/^7/57 



* HFGI TWQI’XCJTG - MR><«HF ( TW9 ) *CJ«G + V A8 #V A 8 /2 , -M I . +X8f ( MR-X 8 I *VR ft WA( 

MIK I*2./( 1. + XP + I mr-X9T-XVR9WA( 

VW<;=VRQWA( IM )*VA9 
A9 = SORT(GC*R*TS<? ) 

AMN9=VA9/A9 

SPECIFIC THRUST AND PROPULSIVE EFFICIENCY CALCULATIONS 

StHRJJS = li t L.TFAR »/ ( l.C+BYD+VAS - Jl.O/Il.OtBYn^QV.* 

*(BY*BY*X9I /( l.+BY)*VA9 - ( BY / ( 1 . +B Y I) *OV *■ ( WGR ( N I -BV*X9 ) / ( 1 . + BY » * 

* VW9 - (WGR(N) l.+BY) )*VWO + STADD) /G 
THRUST=STHRUS*(CWF + FMF) ♦ G 
TP=THRUST«QV 

FTAPRO=TP/(TP«^TMF/2.0*(VA5-OV>#*2 + ( FMF + X9*FMF I /2.* ( VAO-OV) ^*2 «• 

* (HMF-FMFYX9)/2.*( VW9-VW0>**2> 

FUMF=( TMF - CMF » 

FUWF=FUMF*G 

SFC=(F11WF /THRUST) ♦3600.0 
WRITE (6, 20) 

2C FORMAT! ////tT8,» PSR 87* , T18, « VR8WA* .T28t ' TSR8WA* . T38, • VR9WA* , T48, • T 
♦SR9WA' , T6D, • ER » iT69, *WGR' ) _ . 

WRITE(6,21 )PSR87( IJ ) , VR8WA! IK) , T SRP WA (I L ) ♦ VR9W A ( IM) , TSR9WA! IN) .BY, 
*WGR(N) 

21 F0RMAT(2X,7F10.2I 
WP ITEI6.22) 

22 FORMAT (/,Te,'VA8',Tl6, » VW8 • , T2-^ , • VA9* , T32 , ' VW9* ,T4r ,'TS8*,T4B,»T 
♦W8* ,T56, • TS9* ,T64, *TW9* ,T71 , • STHRUS • , T8 1 , • AMNS • , T88 , • AMN9' , T94, • ET 
♦APRO*,T102,'SFC*,T109,'FAR',T116,'X9M 

WRITE(6.23)VA8,VW8,VA9,VW9,TS8,TW8,TS9,TW9,STHRUS, AMN5, AMN9,ETAPP0 
*,SFC,FAR,X9 

23 F0RMAT(4X,8Fe.2,F9.2,5F7.3,F7.4) 

GO TO 26 

24 WRITE(6,25) 

25 FORMAT! T5, 'TURBINE EXIT PRESSURE IS I OWER THAN ATMOSPHERIC) 

26 CONTINUE 

27 CONTINUE 

28 CONTINUE 

29 CONTINUE 

30 CONTINUE 

31 CONTINUE 

32 CONTINUE 

33 CONTINUE 

34 CONTINUE 

35 CONTINUE 

36 CONTINUE 
RETURN 
END 
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